Nonequilibrium phases in hybrid arrays with flux qubits and NV centers 
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We propose a startling hybrid quantum architecture for simulating a localization-delocalization transition. The 
concept is based on an array of superconducting flux qubits which are coupled to a diamond crystal containing 
nitrogen- vacancy (NV) centers. The underlying description is a Jaynes-Cummings-lattice in the strong-coupling 
regime. However, in contrast to well-studied coupled cavity arrays the interaction between lattice sites is me- 
diated here by the qubit rather than by the oscillator degrees of freedom. Nevertheless, we point out that a 
transition between a localized and a delocalized phase occurs in this system as well. We demonstrate the possi- 
bility of monitoring this transition in a non-equilibrium scenario, including decoherence effects. The proposed 
scheme allows the monitoring of localization-delocalization transitions in Jaynes-Cummings-lattices by use of 
currently available experimental technology. Contrary to cavity-coupled lattices, our proposed recourse to styl- 
ized qubit networks facilitates (i) to investigate localization-delocalization transitions in arbitrary dimensions 
and (ii) to tune the inter- site coupling in- situ. 
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I. INTRODUCTION 

In recent years, a variety of novel experimental approaches 
have enabled tests of fundamental quantum physics such as 
superpositions, entanglement, tunneling or quantum phase 
transitions in artificial devices. Prominent examples hereof 
are quantum circuits, quantum dots or optical lattices dill. 
Apart from their fundamental relevance, quantum technolo- 
gies allow realizing quantum information processors that 
bring along the potential of carrying out specific tasks at expo- 
nentially reduced computation time [3]. Furthermore, quan- 
tum simulators of Feynman type (HIH, employed to simulate 
the dynamics of one quantum system by means of another one 
are a pivotal example of quantum speed-up as compared to a 
classical computer. While both systems share the same dy- 
namics, the simulator offers far more configurability and is 
better accessible for a measurement. 

Generally, the possibility of observing quantum effects 
strongly depends on the coherence properties of the under- 
lying system. With regard to superconducting circuit qubits, 
the advantages of versatile manufacturing, detection and ma- 
nipulation are paid for at the price of quite high decoherence 
rates as compared to trapped ions or spin qubits. One way out 
is given by the recently emerging field of hybrid systems (6]- 
fT4ll . The main motivation for building hybrid systems is to 
combine two advantages: the addressability of artificial quan- 
tum circuits and the long coherence times of elemental sys- 
tems such as nitrogen- vacancy (NV) centers in diamond or in 
polar molecules. Usually these hybrid systems are motivated 
by using the natural spins for building quantum memories. 

In this work, we point out an alternative application for 
the exploration of many-body physics such as quantum phase 
transitions. In particular, we investigate the localized and de- 
localized phases that occur in Hubbard-like models such as the 
Jaynes-Cummings (JC) lattice 1 15 - 18 1. In the localized phase, 
excitations are localized at individual lattice sites, whereas 
they are delocalized across the lattice in the delocalized phase. 
We propose an intriguingly simple layout for simulating a JC- 
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Figure 1. (a) Schematic of a single flux qubit coupled to a diamond 
crystal with NV-centers embedded, (b) Energy diagram for the qubit 
with level splitting ujq and a NV-center. By applying an external mag- 
netic field, the level with spin projection m = ±1 becomes resonant 
with the upper qubit level, (c) The JC-array with tunable qubit-qubit 
coupling. Here, adjacent qubits are connected via auxiliary tunable 
coupler qubits. Each of the qubits couples to spatially separated re- 
gions of the crystal. The coupler qubit does not couple to the spins 
because it is far detuned from the qubits and consquently from the 
NV- spins as well. 



lattice: We use the combination of an already experimentally 
well-proven flux qubit array together with a single large NV- 
center crystal. By means of numerical studies, we corroborate 
that this system exhibits localized and delocalized phases. Fi- 
nally, we demonstrate that these phases can be identified by 
monitoring the signatures of the non-equilibrium system dy- 
namics in presence of decoherence and dissipation upon em- 
ploying experimentally accessible parameters. 



II. HYBRID QUBIT-RESONATOR MODEL 

As the elementary unit for the JC-array we propose a hybrid 
combination of a flux qubit and an ensemble of independent 



spins, given here by the NV-centers in a diamond 1 19 , 20 1. As 
illustrated in Fig.[TJa), the spin crystal is placed in proximity 
of a qubit loop. The spin-qubit Zeeman interaction is medi- 
ated by the magnetic field that stems from the qubit 's persis- 
tent currents. A weak external field splits the spin degener- 
acy so as to shift one spin transition into resonance with the 
qubit 's transition frequency ujq. Reducing the spin to its two 
lowest levels with mutual energy spacing ujc (see Fig. [T] (b)) 
and applying the rotating wave approximation (RWA) we can 
express the Hamiltonian of an ensemble of N spins coupled 
to a single flux qubit at its degeneracy point sls (H = 1) 
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(1) 

Here, a^ = a^ ± ia^ and r^ = r^ ± ir^ are raising and low- 
ering operators with respect to the qubit 's (a) and the spins' 
(r/e) Pauli matrices. The coupling strengths gk between the 
qubit and the individual spins are proportional to the magni- 
tude of the qubit's field at the spin positions 1 19, 20 1. Assum- 
ing the case of a large spin ensemble with low polarization, 
i.e., close to its ground state, we introduce a collective opera- 
tor: i.e. with g = {^^ l^feP)^/^ we set a^ = g-^ Y.k 9kr^ 
together with its hermitian conjugate a, yielding approxi- 
mately the bosonic commutation relation, [a^a^] = 1 (See 
Appendix A). Thus, we interpret the ensemble as an effective 
bosonic mode and arrive at the effective Jay nes -Gumming s 
model, reading as 



Hjc = LUqCr^cr 
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It describes a collective harmonic oscillator mode being cou- 
pled to a two-level system with the interaction strength g. The 
collective coupling g is enhanced by a factor of ^/N compared 
to the root mean square of the individual couplings gk, see Ap- 
pendix [A| Recent experiments achieved coupling strengths as 
strong as ^ ^ 27r X 35 MHz Ell. 



III. THE JAYNES-CUMMINGS LATTICE WITH 
QUBIT-QUBIT COUPLING 

A general advantage of superconducting circuits is their 
scalability and the rich variety of coupling mechanisms that 
can be implemented on a chip. In particular, arrays of 
flux qubits with tunable coupling strength between individ- 
ual qubits have been realized using a SQUID or ancilla flux 
qubit ll22H26ll . In recent experiments manipulating coupling 
strengths in situ and the engineering of various types of circuit 
connectivity has become feasible |27 -31 1. In this work, we re- 
strict ourselves to a chain of qubits with tunable nearest neigh- 
bor interaction. This array of qubits can be readily turned into 
a JC-lattice of coupled qubit-oscillator systems by putting one 
NV-center crystal on top, as sketched with Fig. [ijc). As ar- 
gued above, the spin crystal adds an effective harmonic oscil- 
lator degree of freedom to each site of the array. Apart from 
the possibility to tune the coupling between the sites of the 
lattice, the most appealing aspect of this hybrid architecture is 



simplicity. Furthermore, the harmonic oscillators in the form 
of the spin crystal exhibit excellent coherence properties, ho- 
mogeneous transition frequencies and coupling strengths, all 
implemented here within a reduced geometric dimension as 
compared to coplanar waveguide resonators. 

For well-separated qubits, we can neglect their mutual in- 
ductance as well as the cross-coupling of one qubit to the spin 
ensemble of another site, being even one order of magnitude 
smaller. This JC-lattice with M sites is thus described by 
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HjcL = Yl ^^c,j + J ^ «^,+i + h. c.) 



(3) 



with the single-site Hamiltonians Hjcj given in Eq. ([2]). Here, 
J denotes the uniform qubit-qubit coupling strength and the 
operators a^ describing the creation and annihilation of a 
qubit excitation at the j-th site. 

A subtle but salient difference between our model and pre- 
viously studied JC-lattices is in the interaction mechanism be- 
tween individual lattice sites. While we propose an inter-site 
coupling mediated by the qubits, previous works have dealt 
with the complementary approach where the lattice sites in- 
teract via the oscillator degrees of freedom, as in coupled cav- 
ities 115-18, 32, 33] or in superconducting resonators [34, 35|. 
In the latter case, the coupling part of the JC-lattice Hamilto- 
nian ^ assumes the form aja^+i + h.c. 



IV. EQUILIBRIUM PROPERTIES OF THE JC LATTICE 

For Bose-Hubbard-like models the occurrence of a quan- 
tum phase transition between localized and delocalized phases 
has been extensively studied 1361 . Analogous transitions 
have been investigated with polaritons in JC-lattices |15- 
[T8l[32j[34l|35l. Here, the term polariton refers to the eigen- 
states |n, ±) of the single-site Hamiltonian Hjc [Eq. (|2])]. The 
excitation number n, being the eigenvalues of the operator 



A/- 



a^ a 



a~^a , are conserved due to [Hjc^J^] = 0. 



Similarly, the full JC lattice Hamiltonian ^ conserves the 
total number of excitations in the lattice. The ground state 
|0) has no excitations n=0, while the states |n,±) (n>0) 
are each twofold degenerate with respect to A/". If the qubit 
and the resonator are in resonance, uJc = cjOq, the polaritonic 
states are symmetric (+) and antisymmetric (— ) superposi- 
tions |n,±) = (|n-l) It) ±|n)||))/V^ of the oscillator Fock 
states \n) and the qubit ground (||)) and excited (|t)) states, 
respectively. 

The localization-delocalization transition we consider in 
this work takes place for the lowest energy state in the sub- 
space with one average excitation per site |[l7l[T8l. As we ar- 
gue with Appendix[D) for weak inter-site coupling J, no inter- 
conversion between the + and — polaritons occurs l T7ll34ll . 
Therefore, in order to obtain analytical estimates, we can ne- 
glect |n, +) polaritons and restrict our studies to |n, — ) polari- 
tons which are lower in energy. We then introduce the "effec- 
tive repulsion" S = E\2-) — 2£^|i _), i.e.. 
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Figure 2. Transition between the localized and delocalized phases 
in a qubit-coupled JC-array. (a) Fluctuations of the number of ex- 
citations at a certain site var(A/'j) for a two- site setup (both sites 
j = 1,2 yield the same plot) as a function of the inter-site coupling 
J and the qubit-ensemble detuning A. The dark shaded region indi- 
cates that the system in the localized phase, while the brighter areas 
are related to large fluctuations, i.e., the delocalized phase. The two 
side panels depict a horizontal cut along J = O.lg and a vertical cut 
along A = 2 X 10~^g (i.e., very close to qubit-oscillator resonance), 
respectively. There, the solid black curves depict var(A/'j) for the 
two-site setup as in the central panel. For comparison, we have in- 
cluded the fluctuation characteristics var(A/'j) for longer JC-array s 
with A^=3-5 sites in ascending order, where j denotes a central site 
of the array, (b) Comparison of QQ- and CC-coupled chains. In the 
latter, the transition occurs at lower detunings due to higher effective 
coupling of polaritons between adjacent sites. This can be seen in 
the lower plot: Changing the detuning affects both the effective re- 
pulsion S, as well as the effective coupling Jeff. The transition occurs 
when 6 and Jeff cross. 



This positive-valued repulsion increases with the qubit- 
oscillator coupling strength g and decreases with the detuning 
A = ujq — Uc. It measures the extra energy needed to insert 
two polaritons into a single site as compared to distributing 
them across two sites. Thus, a large repulsion promotes an 
even distribution of excitations over the lattice sites. In this 



case, the system eigenstate is approximately given by a prod- 
uct of the local single-site eigenstates |1, — )j. 

By contrast, a large inter- site coupling quantified by J fa- 
vors delocalized excitations, i.e., momentum eigenstates that 
are given by a superposition of product states, each with dif- 
ferent n. Thus, modifying J or the repulsion (e.g. by means 
of A) one ends up in two extreme regimes: the localized or 
the delocalized phase. 

As stated at the beginning of this section, the Hamilto- 
nian ([3]) conserves the total number of excitations. Therefore, 
the fluctuation of the excitation number in a particular lattice 
site, var(A/'j) = (J^f) — (ATj)^, is used as an order parame- 
ter in JC-lattices I.34J . For var(A/'j) = 0, the excitations are 
trapped, and the system is in a localized phase. By contrast, 
large fluctuations indicate the delocalized phase. 

In order to investigate the transition between the localized 
and the delocalized phases numerically we calculate var(A/'j) 
for a setup with two sites. As indicated in the main panel of 
Fig. [2] (a), a transition between the two different regimes char- 
acterized by zero and finite var(A/'j) occurs upon a change of 
A or J. This two- site setup already exhibits the same qual- 
itative features as longer arrays of finite length, as we cor- 
roborate in the side panels of Fig. [2ja). There, we compare 
the variance for arrays with two or more sites by means of 
two cross sections through the main panel for fixed values of 
A and J, respectively. Thus, the elementary two-site setup 
- readily feasible with present-day experimental techniques - 
already allows for a good qualitative estimate of the transi- 
tion properties of a JC-array. Experimental feasibility will be 
further discussed in section Ivll 

To gain analytical insight, we express the coupling be- 
tween the individual sites in terms of the relevant polari- 
tonic basis states |n, —)j. In doing so, we can approximate 
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. |n + 1, — )j(n, — Ij, where the coefficients 



Sn, depend on A, ^ and n and their explicit form is detailed 

in Appendix [c] Thus, two sites initially in the state |1, — ) are 
coupled with the effective strength (see Appendix [C]) 



-^eff = -^-^O, Si 



(5) 



In the lower panel of Fig.[2tb), we compare Jeff to the effective 
repulsion strength S [Eq. Q], both plotted as functions of A at 
fixed J. We find that the observed crossing point of Jeff and S 
closely matches the location of the localization-delocalization 
transition. 

Furthermore, we compare our results for a JC-array with 
qubit-qubit (QQ)-coupling to a similar setup with cavity- 
cavity (CC)-coupling, i.e., an array in which the individ- 
ual sites interact via their oscillator degrees of freedom 
J(aj a j+i +h.c). For this latter scenario, we find that the tran- 
sition to the delocalized phase already occurs at smaller A, see 
Fig.[2jb), top panel. As in the (QQ)-coupled case we calculate 
Jeff via the relevant polaritonic basis states and indeed find a 
larger effective inter- site interaction that hence explains the 
observed transition point. The interested reader can check the 
explicit coupling coefficients in Appendix |C] With increas- 
ing detuning, the |n, —)j polaritons become more and more 
bosonic, i.e. only the oscillator degree of freedom is excited, 
|n, —)j ^ \i)j\'^)j- This allows for a simple explanation of 



the different trends of Jeff in Fig.[2jb) when increasing the de- 
tuning. In a qubit-coupled array, the bosonic excitations must 
hop via the route oscillator-qubit-qubit-oscillator to reach the 
next lattice site, therefore Jeff is small. By contrast, in the 
cavity-coupled setup the bosonic excitations can hop directly 
to the next-site oscillator and Jeff is therefore larger when the 
excitations are purely bosonic. 



V. NON-EQUILIBRIUM DYNAMICS 

The characterization of localized and delocalized phases at 
equilibrium is helpful in exploring the physics in JC-lattices 
with qubit-qubit coupling and contrasting it with JC-lattices 
that interact with a cavity-cavity coupling. However, the as- 
sumption of staying in the subspace with one mean excitation 
per site is not completely realistic in practice. In particular, 
non-equilibrium processes such as dissipation and decoher- 
ence are crucial in solid state devices. Next we point out that 
the signatures of the localization-delocalization transition re- 
main preserved even in the presence of dissipation. A cor- 
responding measurement only requires state preparation and 
qubit readout. 

We model dissipation for both the qubit and spin ensemble 
by means of a quantum master equation, which for the JC 
lattice assumes the form (at zero temperature) |[37l[38l 



m 



[HjcL. q]^Y1 (^^^«. [^(^)] + ^^^-- [^(^)]) 



(6) 

The Lindblad dissipators Lq act on the density operator p as 
Lo[p] = OpO^ - \{0^0p + qO^O). The operators O = 
{cr ~, ttj} describe the system-bath coupling of the j-th qubit 
and oscillator, respectively, while 7^ and 7c are the associated, 
uniform decoherence rates. 

The system is initially prepared with one |1,— ) polari- 
ton in each site. Calculating the system dynamics numer- 
ically, we obtain the time-dependent probability P2{'t) = 
Tr{n2^(t)} of finding two excitations in one site, where 
n2 = |2)|;)(^|(2| + |l)|t)(e|(l|. If the system is in the local- 
ized phase we expect P2{t) to remain close to zero. By con- 
trast, in the the delocalized phase, P2 {t) reaches finite values 
over time. This behavior is depicted in Fig. [Sj^b). Here, it is 
also visible that the system evolves eventually into its ground 
state due to decoherence. 

In order to quantify the phase we introduce the averaged 
probability 



dtP2{t). 



(7) 



In order to take into account the dynamics before relaxation 
into the ground state dominates, the integration time should 
< 



fulfill T > J-^^ but T < min{7-\7-i}. Fig.jsja) depicts 
P2 as a function of both the hopping parameter J and the de- 
tuning A similar to the equilibrium analysis. For compari- 
son, the white dashed line marks the parameter regime where 
the phase transition occurs in the equilibrium case in Fig. [2] 



0.15 




0.05 



^0.2 

^0.1 













/ \ 


— A = O.lg, J = 0.1g- 


; \ /x 


-- A =.O.lg, J =lg _ 


/ 1 / » /N 


_ 


/ '.' \/ \.'~V .'>. 





2 3 

t (units of 7^"*^) 



Figure 3. Non-equilibrium signature of the phase transition, (a) 
Time- averaged probability A to find two excitations in a single site. 
The coupling is assumed as 5^ = 27r x 10 MHz and the decay rates of 
qubits and oscillators are 7q = 27r x 1 MHz and 7c = 27r x 0. 1 MHz. 
We choose the integration time as T = 67"^. The dotted line marks 
the boundary (1/2 max{var(A/'j)}) where the phase transition occurs 
in the equilibrium case of Fig. [2] (b) Time evolution of P2 in two 
exemplary points in the delocalized (dashed line and dashed cross in 
(a)) and localized phase (solid line and solid cross in (a)), respec- 
tively. 



While we find a good agreement for small values of A, the 
border between both phases is not resolved in the far-detuned 
limit where the effective inter-site coupling strength Jeff de- 
creases below the decoherence rates {7^, 7c}. 



VI. EXPERIMENTAL FEASIBILITY 

We next address the feasibility of the JC array proposed 
in this paper. The case of strong coupling between an en- 
semble of NV-centers (ujc ^ 27r x 2.88 MHz) and a flux 
qubit has previously been reported experimentally with cou- 
pling strengths up to ^ ::^ 27r x 35 MHz II2TII . Thus, using 
our estimations, we can safely consider a coupling strength of 
g ^ 27r X 10 MHz. On the other hand, the experimental acces- 
sible tunable qubit-qubit couplings are between 27r x 1 MHz < 
J < 27r X 100 MHz ||23|[25||26l. AH together sets the oper- 
ation range to 0.1 < J/g < 10. We next consider the decay 
rates of the involved subsystems. Realistic values for the qubit 
decay rates are 7^ = 27r x 1 MHz (i.e. jq/g = 0.1) and for the 
spin decay rates are 7c < 27r x 0.1 MHz (i.e. jc/g < 0.01). 

These values were used to obtain the plots in Fig. [3] There- 
fore, current technology allows for monitoring both phases. 



Finally, readout of the number of excitations at a specific site 
can be performed by measuring the qubit dynamics. 



VII. CONCLUSIONS 

In summary, we have introduced a novel JC-lattice based 
on a hybrid combination of flux-qubits and NV-centers. In 
contrast to JC-lattices based on coupled cavities or supercon- 
ducting resonators, the harmonic oscillator degree of free- 
dom ("cavity") is smaller in size than the qubit ("atom"). 
This allows one to couple the individual JC- sites via the 
qubits instead of the harmonic oscillators. We have argued 
that similarly to cavity-coupled JC-lattices a localization- 
delocalization transition can be observed in these novel qubit- 
coupled JC-lattices. Even though localization-delocalization 
transitions in JC-lattices have been proposed theoretically 
some time ago, they could not be observed in an experiment 
yet. Our proposal relies on a straightforward modification of 
already realized flux qubit arrays by simply mounting a sin- 
gle NV-center crystal on top. This minimal modification of a 
common setup opens the possibility of studying many-body 
phenomena in strongly coupled hybrid architectures within 
state-of-the-art experimental technology. Apart from its sim- 
plicity, further advantages are the possibility to investigate 
localization-delocalization transitions in arbitrary (even frac- 
tal) dimensions and to tune the inter- site coupling in- situ by 
using common techniques for building flux qubit networks. 
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Appendix A: Collective modes 

Here we outline the approximation that allows us to express 
the spin ensemble by a collective bosonic operator. With ho- 
mogeneous coupling for all the spins one can express the spins 
by a collective angular momentum operator and then apply a 
Holstein-Primakoff approximation to yield a bosonic opera- 
tor. However, for inhomogeneous couplings this is not readily 
possible because the collective operator does not fulfill angu- 
lar momentum algebra. Nevertheless, we can arrive at col- 
lective bosonic operators: We start with the Hamiltonian of 
spins (inhomogeneously) coupled to a qubit with the individ- 
ual coupling strength g^ 
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H = ujqa^a ^uJc^TjItj^ ^^{gkcr r+ + h. c.) . 

k k 

(Al) 



Here, r^ are the Pauli raising and lowering operators of the 
spins in the ensemble and a^ the ones of the qubit. The spins 
are taken to have the homogeneous energy splitting ujc and 
the qubit the splitting ujq. We next introduce the collective 
operator 
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where g is the root mean square of the individual couplings, 
g'^ = ^j^\gk\ /N. In the low polarization limit, where al- 
most all spins are in the ground state, it follows that these op- 
erators approximately fulfill bosonic commutation relations. 
To validate this we calculate the commutator 

^ ki 

1 v-^, 



N-g-^ 



k|'(lfc-2T+T- 



Inserting the definition of g yields 

2 
Wg 



[a, at] =1-^^15,1 V+T,-- (A3) 



For states with only a few spins excited the second term is <C 1 
and thus a^ and a obey approximately bosonic commutation 
relations. 



[a, a^] ^ 1 . 



(A4) 



Using these collective operators the coupling term of Hamil- 
tonian \A\) becomes 



(A5) 



7V^(cr+a + cr"a'^) 



If we start with the spin ensemble in the ground state |0) and 
with the qubit excited, Rabi oscillations can transform exci- 
tations into the ensemble that assume the form of generalized 
Dicke states 
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5fc|0i...lfe...0jv) = |l). (A6) 



Here |0i . . . 1^ . . . Oat) denotes a state where all spins are in 
the ground state except the k-\h spin. Higher excited states 
are defined by 



^(«Tlo). 
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If we restrict our Hilbert space to the set of states with this 
symmetry, we can express ^^ ^t^k ^^ ^^^ ^^^^ P^^ ^^ 



Hamiltonian ( |A1| ) by collective operators as well. To this end, 
we show that ^^ ^k'^k gi^^s the number of collective exci- 
tations for state |0), 



The individual Jaynes-Cummings energies E^^a at a certain 
site are thus given by 
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^r+r,-|0)=0|0) 



(A8) 



and we use induction to show that if it is trae for \n) it is also 
true for |n+ 1): 
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En,± =nuJc^-± Vng^^Ay^ , (B2) 

With the detuning A = Uq — Uc. Note that the |n, — ) po- 
laritons are lower in energy than their |n, +) counterpart. The 
eigenstates are 

|n, -) = cos6>n|n,|) -smOn\n- l,t) 

|n,+) =sinl9n|n,;)+cos6>n|n- l,t) (B3) 

and |0) with Eq = where the mixing angle On is defined as 

1 



On = - arctan I , , 
2 V^/2 



(B4) 



Vn 



= n\n + 1) + l|n + 1) = (n + 1) |n + 1) . 

In the next to last step in Eq. ( |A9| ) we insert the relation: 

S/c ^k^k\^) ~ ^1^)' valid by induction. Besides, we 
use the Pauli matrices conmutation relation: [t^,t^^] = 
^ki (ife — ^^^k^k) ^^^ notice that Tj^Tj^ = 0. Therefore, we 
can write in the subspace of the collective excitations | n) 



With increasing detuning the |n, — ) polaritons resemble more 
and more pure bosonic excitations, |n, — ) ^ |^) |i), while the 
|n, +) polaritons exhibit an excitation in the qubit, |n, +) ~ 

|n-l)|t>. 



Appendix C: Hopping term in the polaritonic basis 

Next we need to express the inter-site-hopping terms. 
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In conclusion, restricting ourselves to the Hilbert space of 
the states \n) the initial Hamiltonian ( |A1| ) can be recast as 



in the local JC-basis basis as well. The operator aj~ takes the 
form 

CX) 

cr + = ^ Y^ Sna(3\n^l,f3)j{n,a\j , (C2) 

n^0a(3=± 



H = w^cT+a- + w,a^a + ^v^ {a+a + cj- a^) . (Al 1) ^^ere the coefficients s„„^ are given by 



Besides, in the low polarization limit a^ and a obey approx- 
imately bosonic commutation relations. Finially, It follows 
that the collective coupling is enhanced by a factor of \fN 
as compared to the root mean square of the couplings to the 
individual spins. 



Appendix B: The polaritonic basis 



Sn Sn-+ \ ^ / -COSe/nSm6/n+l COS f^n COS t/n+l 

Sn+- 5n++ ) ~ \ -sin6>nsin6>n+i sin l9n cos 6>n+l 

(C3) 
These coefficients depend on the coupling g and the detun- 



ing A between the oscillator and the qubit via ( |B4[ ). Note 
that the operator a^ acts on a polariton state |n, a = ±)j by 
transforming it to a linear combination of polaritons ± with 
an additional excitation 



The eigenstates of a single JC-Hamiltonian \A\\\ are the 
so-called polaritons and are denoted by |n, ±)j. Here, n 
describes the number of excitations at a site (sum of qubit 
and bosonic excitations), while the sign defines the polariton 
"species". We now use these polaritons to express the Hamil- 
tonian of the whole JC-chain. 

Using the polariton eigenstates |n, ±) j at site j the individ- 
ual uncoupled JC-Hamiltonians are diagonal. 



(j^\n,a)j = Y^ Sna(3\n^l,i 

/3=± 



/J ' 



(C4) 



H: 



JC 



2^ 2_^^n,a\n,a)j{n,a[ 



(Bl) 



n—O a— da 



The hopping term can now be rewritten in the polaritonic basis 
as 

i7]^°P = J^^^8na/3|n + l,/3),(n,a|, X (C5) 

nn' OLOi' 1313' 

Sn'cx'^' \n', a)j+i {n + 1, /3'| j+i + h. c. 

The hopping term therefore allows transitions between polari- 
tons of different types on adjacent sites. 



1. Effective coupling Je// in the case of two sites 

For a two site array initially in the state |11) = |1,— )i <S) 
|1, — )2 we can derive an effective coupling constant between 
the two sites. Discarding transitions between different polari- 
ton species the dynamics are spanned by {|11), |22)} with 

|22) = I/V2 (|2, -)i |0, -)2 + |2, -)i |0, -)2) . 

(C6) 
In this case the hopping term simply reduces to 

(ll|i7]^"P|22) = J(8i__5o— +5i__8o— ) . (C7) 

Therefore the effective coupling emerges as Jeff = 
J Si So 



2. Effective cavity-cavity coupling 

As in the qubit-coupled case we rewrite the harmonic oscil- 
lator operators in the new basis 



^ tnai3\n^l,f3)j{n,a\j , (C8) 

where the coefficients tna(3 are given by 



E 

n=0a(3^± 



tn = COS On cos ^n+1 a/ti + 1 + sin On sin On+1 ^/u 

tn \- = cos On sin On+w/n + 1 — sin On cos ^n+1 V^ 

tnH = sin On cos On-\-iVn + 1 — cos 6>n sin ^n+1 V^ 

tn-\-\- = sin 6>n sin On-\-i Vn + 1 + cos 6^^ cos ^^+1 V^ 



For a two site cavity coupled array the effective coupling is 
finally given by Jg// = Jh to 



Appendix D: Conservation of polariton type 

In the previous section we have seen that, generally, the 
coupling transfers polaritons from one site to another and may 
change the polariton type (±) on the sites involved. 

However, not all of these possible transitions actually have 
to occur. Transforming the coupling terms to the interaction 
picture with respect to the uncoupled Hamiltonian; i.e.. 



Ho = YHf = V En,a\n,a)j{n,a[ 



j,n,a 



(Dl) 



we detect that the individual terms of the coupling Hamilto- 
nian, Eq. ( [C5| ), acquire time-dependent rotating phases 



4^nn'aa'(3(3'{^) 



(D2) 

We note that the overall frequency of the time dependent 
phase - given in terms of the energies of the states it is cou- 
pled - broadly varies in magnitude. Thus, under the premise 
of a small inter-site coupling J, we can identify fast oscillat- 
ing terms and neglect those within a rotating wave approx- 
imation (RWA). Particularly, terms connecting states which 
are far apart in energy (:^ J) and thus are rotating fast, are 
the ones that can be neglected. 

This treatment is certainly useful in the case of large detun- 
ing or strong coupling strength g. Then, polaritons of differ- 
ent sign lie far apart in energy. Thus, their interconversion can 
safely be neglected and the initial polariton species (±) is thus 
conserved. 
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